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Electron transport in mesoscopic conductors has traditionally involved investigations of the mean 
current and the fluctuations of the current. A complementary view on charge transport is provided 
by the distribution of waiting times between charge carriers, but a proper theoretical framework for 
coherent electronic systems has so far been lacking. Here we develop a quantum theory of electron 
waiting times in mesoscopic conductors expressed by a compact determinant formula. We illustrate 
our methodology by calculating the waiting time distribution for a quantum point contact and find a 
cross-over from Wigner-Dyson statistics at full transmission to Poisson statistics close to pinch-off. 
Even when the low-frequency transport is noiseless, the electrons are not equally spaced in time due 
to their inherent wave nature. We discuss the implications for renewal theory in mesoscopic systems 
and point out several analogies with energy level statistics and random matrix theory. 

PACS numbers: 72.70.+m, 73.23.-b, 73.63.-b 



Introduction. — The distribution of waiting times be- 
tween elementary physical events is an important tool 
to characterize and investigate temporal correlations and 
transport statistics of stochastic processes [1]. Waiting 
time distributions (WTDs) play a significant role in var- 
ious branches of science and technology, for instance in 
quantum optics [2 , 3 and single- molecule chemistry [4]. 
There, the dynamics is often Markovian and renewal the- 
ory applies such that the WTD fully defines a random 
walk in time and most statistical properties can be ob- 
tained from the WTD alone [5]. In contrast, for systems 
where renewal theory is not valid, the situation is more 
complex and various statistical quantities provide differ- 
ent information about a stochastic process. The dynam- 
ics is then similar to a random walk with memory effects. 
Investigations of WTDs have until now mostly dealt with 
either classical systems or Markovian quantum systems. 

In electron transport, the waiting time between charge 
carriers has only received limited attention. One ex- 
ception concerns weakly coupled nano-structures whose 
dynamics essentially is classical [Gj |7j . In the opposite 
regime of fully coherent conductors, a theoretical frame- 
work for WTDs has so far been missing. The attention 
has instead been devoted to the full counting statistics 
(FCS) of transferred charges which for non-interacting 
electrons is elegantly expressed by the Levitov-Lesovik 
determinant formula [8j. Typically, FCS concerns the 
limit of long times and important short-time physics may 
be lost (see however Ref. [9 and references therein). To 
describe short-time physics and correlations, WTDs are 
particularly useful and a theory of WTDs in mesoscopic 
conductors is thus of fundamental importance. More- 
over, with the rapid progress in single-electron detection, 
there is hope that measurements of WTDs for a meso- 
scopic device could soon be within reach. 

In this Letter we develop a quantum theory of WTDs 
for mesoscopic conductors. The generic non-equilibrium 




FIG. 1. (color online). Quantum point contact (QPC) con- 
nected to voltage-biased electrodes. The QPC transmission 
is denoted as T, the Fermi energy is Ep^ and V is the voltage 
applied to the source electrode. The average time separation 
of the in-coming electrons is f = h/eV. The distribution of 
waiting times r between transmitted electrons is determined 
by the many-body state of the in-coming electrons as well as 
the QPC which may cause electrons to reflect back. Reflected 
(missing) electrons are indicated by dashed lines. 



system consists of electronic leads connected to a scat- 
terer described by the scattering matrix S. As a cen- 
tral result we formulate a compact determinant formula 
for the WTD in terms of the transmission amplitudes in 
S. We employ a first- quantized many-body description 
which explicitly incorporates the fermionic statistics of 
the in-coming particles imposed by the Pauli exclusion 
principle [10]. This is a generally applicable methodol- 
ogy. We illustrate it by calculating the WTD of a voltage- 
biased quantum point contact (QPC), Fig. [l] For a fully 
open QPC, the low- frequency transport is known to be 
deterministic and regular [11 . Nevertheless, we find that 
the electrons are not equally spaced in time due to their 
inherent wave nature. Moreover, unlike the FCS, which is 
always generalized binomial [TO, 12^ , the WTD exhibits 
a cross-over from Wigner-Dyson statistics for an open 
QPC to Poisson statistics close to pinch-off. These find- 
ings show how the WTD and the FCS provide separate, 
yet complementary, views on quantum transport in meso- 
scopic conductors. In general, renewal theory does not 
apply for mesoscopic conductors and no simple relations 
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exist between the WTD and the long-time FCS. We con- 
clude by pointing out a number of appealing analogies 
between WTDs and energy level statistics in random ma- 
trix theory J4]. 

Formalism. — The distribution of waiting times r be- 
tween charge carriers is denoted as W(r). It can be ex- 
pressed in terms of the idle time probability n(r): Given 
a random time to, the probability that no electrons are 
detected in the time interval [to, to + r] with r > is 
denoted as n(r). For stationary processes n(r) does not 
depend on to- The WTD can be written 

W{r) = {T)^niT). (1) 

Here (r) is the mean waiting time which ensures the 
proper normalization drW{r) = 1, since 11(0) = 1 
by definition and n(r) decays to zero at large times. 
This important relation can be derived by picking a ran- 
dom time to and assuming that the last prior detec- 
tion of an electron occurred at the earlier time tg < to- 
The idle time probability is then n(r) = J^^dte[l — 

/^*°+^ dsW{s-te)]/{r), where the square brackets contain 
the probability that no electrons are detected in the time 
interval [te,to + r]. We also integrate over all possible 
times for the last electron to be detected, — oo < tg < to, 
using that detection events are uniformly distributed in 
time with weight 1/(t). Finally, a suitable change of in- 
tegration variables yields n(r) = du dvyV{v)/ (r) 
which explicitly shows that the idle time probability is 
independent of to and immediately leads us to Eq. Q. 
Importantly, in deriving Eq. ([T]) we have made no Markov 
assumption nor relied on any renewal property. 

Our next task is to evaluate the idle time probabil- 
ity for a generic mesoscopic conductor described by the 
scattering matrix S. To this end we make use of the 
first-quantized many-body formalism developed by Has- 
sler et al. [10] . Our methodology applies to arbitrary ge- 
ometries of the electronic leads and the scatterer, but to 
keep the discussion simple we consider a one-dimensional 
problem with a single in-coming and out-going quan- 
tum channel connected to the scatterer. The electrons 
are non-interacting and spinless and the temperature is 
zero. An applied voltage bias V brings the system out 
of equilibrium such that in-coming electrons on the left 
side of the scatterer occupy states in the energy window 
ly = [Ef, Ep^eV] above the Fermi level Ep. The scat- 
tering states take the form (^k{x) = e*^^ +r/ee~*^^, x < 0, 
and (^/e(x) = t^e*^^, x > > 0, where Vk and tk are the 
refiection and transmission amplitudes in S and k is the 
momentum. We need not specify the scattering states in 
the scattering region [O,^^]. 

We work close to the Fermi level, where the dispersion 
relation E^ = hvpk is linear in k and all wave components 
propagate with the Fermi velocity vp. Next, we split the 
energy window ly into N intervals of size eV/N ^ where 

is the number of incoming particles, and associate to 



the m'th energy interval a time-dependent wave func- 
tion (t)m{x,t) = {x\(l)m{t)) = Jj^dke-'''^^^(pk{x)/V2^. 
Here Im = /^[m — 1, m] is the integration interval of width 
hz = eV/NhvF' All single-particle states are filled and the 
A^-particle many-body state is given by the Slater deter- 
minant l^y^^(t)) of the states \(j)m{t))^ m = 1, . . . , A". 

The detection of individual particles to the right of the 
scatterer, at xq > x^, is effected by the single-particle 
projection operator Q^. Its expectation value with re- 
spect to a single-particle state is the probability of de- 
tecting the given particle in the time interval [to, to + r] 
or equivalently, given the linear dispersion, to detect the 
particle in the spatial range I^- = [xq^xq + vpr]. We 
can therefore write Qr = Jj dx\x){x\ [10 . This is a pro- 
jective measurement. The idle time probability for the 
Slater determinant is then U{r) = (^^^V)! - 
QT-)l^y^^(r)). The expectation value of such a product 
of single-particle operators with respect to a Slater de- 
terminant can itself be written as a determinant and we 
thereby obtain 

H(r) = det(l-Q,). (2) 

The single-particle matrix elements of Q^- are [Qr]m,n = 
{^m{r)\Qr\Mr)) = Jj^dk' Jj^dktltkKrik - k')/27rK 
with the kernel Kr{q) = 26-''^''^^^'^ sm{qvFr /2) /q. Fi- 
nally, to obtain a stationary process we take the limit 
A" ^ oo for which ^ 0. We then have [Qr]m,n — 
hzf^^^t^riKrii^i'^ — m))/27r for m,n = 1,...,A^. Equa- 
tions ([l][2| constitute the central results of this section 
as they allow us to calculate the WTD for an arbitrary 
voltage-biased scatterer. 

Quantum point contact. — We illustrate our method 
by calculating the WTD for a QPC with energy inde- 
pendent transmission |t^rnP = T. Figure [2^ displays 
WTDs calculated for different transmission probabilities 
T. For a fully open QPC, T = 1, the WTD refiects 
only the fermionic statistics and correlations in the many- 
body state. In particular, W(0) = 0, since two elec- 
trons cannot occupy the same state. This suppression is 
similar to the Fermi-hole in the density-density correla- 
tion function of a free electron gas [15 . The fermionic 
correlations also force the WTD to decay strongly be- 
yond a few mean waiting times, where it essentially van- 
ishes. The mean waiting time is (r) = 1/Tt{Qt-=o) with 
=drQr- It easily follows that Tr(Q^^o) = Tr(Q^) = 
{eV/h) X)m=i l^/^mP/A^ = {I)/e is the average (particle) 
current. For the QPC, we find (r) = f/T with f = h/eV. 
This is in line with the common wave-packet picture that 
electrons emitted from a const ant- volt age source into a 
quantum channel on average are separated by the fun- 
damental time scale h/eV [13 . However, the electrons 
are not equally spaced in time due to their inherent wave 
nature which introduces fiuctuations in the waiting time. 
Interestingly, for T = 1 the WTD is well approximated 
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FIG. 2. (color online). Electron waiting times for a quantum point contact (QPC). a) Waiting time distributions (WTDs) 
for different transmission probabilities T of the QPC. The waiting time r is given in units of f = h/eV or (r) — f IT (inset). 
For a fully open QPC, T = 1, the WTD is well approximated by a Wigner-Dyson distribution, Eq. ([3|. Close to pinch-off, 
T :^ 0, the WTD approaches Poisson statistics given by the exponential distribution W(t) ^ e~^/^^^/(T). b) WTD in the 
low-transmission regime, T ^ 0.1. The short-time behavior given by Eq. Q is indicated with a dashed line. Thin lines show 
expansions of the WTD to first, second, and third order in T. Arrows indicate small oscillations with period f. 



by the Wigner-Dyson distribution [14] 

WwD(r) = ^e-4-V-^^ f= A. (3) 

This can be understood using analogies with energy level 
statistics and random matrix theory as we discuss below. 

As the transmission T is reduced below unity, the 
WTD widens as electrons may now reflect back on the 
QPC, allowing for longer waiting times between transmit- 
ted charges. Fig. |2^. Close to pinch-off, T 0, the mean 
waiting time becomes much longer than f and the WTD 
approaches an exponential distribution corresponding to 
Poisson statistics involving rare and uncorrelated tun- 
neling events. In this regime, the WTD is mainly deter- 
mined by the scatterer, although the complete suppres- 
sion at r = persists. In fact, for arbitrary transmission, 
we can expand the WTD to lowest order in r as 

W(r)^y-(-) , r«f. (4) 

The Pauli exclusion principle manifests itself in a r^- 
dependence of the WTD for short waiting times and a 
complete suppression at r = 0. Equation Q is indicated 
in Fig. (2)3 and is in good agreement with the WTD for 
arbitrary transmissions T as we have checked. 

The WTD exhibits a cross-over from Wigner-Dyson 
statistics to Poisson statistics with decreasing transmis- 
sion. The intermediate regime displays small oscillations 
with period f due to the quasi-regular train of incoming 
electrons superimposed on an exponential decay, Fig.|2]3. 
Physically, the wave packets have a large overlap leading 
to small-amplitude oscillations. This situation resembles 
the density-density correlations of a liquid whose parti- 
cles have a large degree of freedom and where the struc- 
ture of the many-body state decays fast with the inter- 
particle distance. In contrast, for a solid-like system. 



whose constituents are placed in a nearly perfect "crys- 
tal" , long-range correlations lead to well-isolated periodic 
peaks in the WTD. This behavior occurs in periodically 
driven classical systems [7 and is also expected for nearly 
non-overlapping quantum states produced by coherent 
single-electron sources [16]. 

To further understand the low transmission regime, we 
expand the idle time probability in T. To this end we 
rewrite it as n(r) = det(l-Q^) = exp [Tr(log{l - Qr})]- 
Since Q^- is proportional to T, we can expand the idle 
time probability as n(r) = 1 - Tr(Q^) + [Tr(Q^)2 - 
Tr(Q^)]/2 + (9(T^). The second derivative with respect 
to time is readily performed, allowing us to express the 
WTD as a series in T. Figure (2)3 illustrates how the 
series improves as more terms are included. In partic- 
ular, to fully account for the long-time behavior of the 
WTD high-order terms must be included. To lowest or- 
der in T we find W(r) = {I)g^'^\T)/e + (^(T^), where 
g'^'^\r) = 1 — s\v? (nt / f) / (lit / fY is the two-point cor- 
relation function [Mj. This is consistent with renewal 
theory [3 . High-order terms, however, cannot be ob- 
tained from renewal theory. The break-down of renewal 
theory becomes evident by considering the cumulants of 
the WTD. 

Moments & cumulants. — The first three cumulants 
of the WTD are the mean ((r)) = (r), the variance 
((r^)) = ((r - (r))2), and the skewness ((r^)) = ((r - 
(r))^). The cumulants of the FCS are similarly denoted 
as ((n^)), m = 1,2,3,..., where n is the number of 
transferred charges. The corresponding Fano factors are 
= ii'^^)) / ii'^)) ^ which for a QPC are well-known and 
readFs = l-TandFa = (1-T)(1-2T) [8^. Importantly, 
for systems where renewal theory applies, the FCS can 
be directly related to the WTD and the Fano factors read 
Fi = {{r'))/{{rr and FJ = 3((r2))V{(rr - {{r'))/{{Tr 



4 




FIG. 3. (color online). Fano factors as functions of QPC 
transmission T. The exact results F2 — 1 — T and F3 = 
(1 — T)(l — 2T) are shown together with predictions based on 
renewal theory = ((r'))/((r))' and = 3((t')) V((t))^ - 
((t^))/((t))^, where ((r^)) are the cumulants of the waiting 
time T. Renewal theory only holds in the low-transmission 
regime (shaded area), where Fm — F^. 



in terms of the waiting time r [3[T7] • These relations pro- 
vide us with a direct test of renewal theory. For a fully 
open QPC, the second Fano factor is zero, correspond- 
ing to regular transport without zero-frequency fluctua- 
tions. However, the width (or variance) of the WTD is 
finite. Fig. |2^, explicitly demonstrating the break-down 
of renewal theory due to the fermionic statistics of the 
in-coming electrons. In Fig. |3] we examine the validity 
of renewal theory as a function of the QPC transmis- 
sion. Only close to pinch-off (T 0), where the trans- 
port is poissonian, we find Fm — F^ in agreement with 
renewal theory. In general, we expect that renewal the- 
ory for mesoscopic conductors is only valid in the low- 
transmission regime where tunneling events are rare and 
nearly uncorrelated. 

Random matrix theory. — There are several analogies 
between WTDs and level spacing distributions in spectral 
statistics ^ . The latter are useful to discriminate com- 
plex systems according to their symmetry class or under- 
lying classical dynamics and they may carry signatures 
of quantum chaos according to the Bohigas-Giannoni- 
Schmidt conjecture For a fully open QPC the anal- 
ogy is due to the equivalence between the ground state 
of one-dimensional fermions and the joint probability dis- 
tribution of eigenvalues of random matrices [TH |19j [20] . 
In particular, the canonical ensemble of random ma- 
trices, labeled by their symmetry parameter /3, can be 
mapped onto the Calogero-Sutherland model of inter- 
acting one-dimensional fermions with coupling constant 
proportional to /3 — 2 Free fermions correspond to 
/3 = 2 if we replace the N coordinates by the N eigen- 
values of a random matrix from the gaussian unitary en- 
semble. As a direct consequence, all spatial correlation 
functions are given by the energy correlation functions in 
random matrix theory. Additionally, since the dispersion 
relation is linear, the system is invariant under Galilean 
transformation and the spatial correlations are identical 



to the temporal correlations. Finally, for an ensemble of 
random hermitian matrices of rank with A" ^ 00, the 
Wigner-Dyson surmise is a good approximation of the 
level spacing distribution and it also agrees well with the 
WTD for a fully open QPC as we found in Fig. [2^. 

Conclusions. — We have presented a quantum theory 
of electron waiting time distributions (WTDs) in meso- 
scopic conductors expressed by a compact determinant 
formula. The WTD of a quantum point contact (QPC) 
exhibits a cross-over from Wigner-Dyson statistics for a 
fully open QPC to Poisson statistics close to pinch-off. 
We explicitly demonstrated the break-down of renewal 
theory for mesoscopic conductors. Finally, we discussed 
analogies between WTDs and level spacing distributions 
in spectral statistics. Open questions to address in future 
work concern the influence of spin and finite-temperature 
effects as well as the WTDs of more complex scatterers. 
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